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ABSTRACT 

XMM-Newton observations of 29 high rcdshift (z>2) quasars, including seven 
radio-quiet, 16 radio-loud and six Broad Absorption Line (BAL) objects, are presented; 
due to the high redshifts, the rest-frame energy bands extend up to ^ 30-70 keV. Over 
2-10 keV, the quasars can be well fitted in each case by a simple power-law, with no 
strong evidence for iron emission lines. The lack of iron lines is in agreement both with 
dilution by the radio jet emission (for the radio-loud quasars) and the X-ray Baldwin 
effect. No Compton reflection humps at higher energies (i.e., above 10 keV in the rest 
frame) are detected either. Over the broad-band (0.3-10 keV), approximately half 
(nine out of 16) of the radio-loud quasars are intrinsically absorbed, with the values 
of Nh generally being 1-2 x 10^^ cm~^ in the rest frames of the objects. None of the 
seven radio-quiet objects shows excess absorption, while four of the six BAL quasars 
are absorbed. The radio-loud quasars have flatter continuum slopes than their radio- 
quiet counterparts (Frl ^ 1.55; Frq ~ 1.98 over 2-10 keV), while, after modelling the 
absorption, the underlying photon index for the six BAL quasars is formally consistent 
with the non-BAL radio-quiet objects. 
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1 INTRODUCTION 

Approximately 10 per cent of all active galactic nuclei 
(AGN) are classified as being radio-loud; that is, the log- 
arithm of the ratio of the radio (5 GHz) to optical (B-band, 
4400 A) fluxes, Rl, is > 1 (Wilkes & Elvis 1987; Kellerman 
et al. 1989). In general, for a given optical luminosity, the X- 
ray emission from radio-loud quasars (RLQs) is about three 
times greater than that from their radio-quiet counterparts 
(RQQs; Zamorani et al. 1981; WorraU et al. 1987), which 
allows them to be studied out to higher redshifts. The high 
effective area of the EPIC (European Photon Imaging Cam- 
era; Turner et al. 2001; Striider et al. 2001) instruments 
on-board XMM-Newton is enabling the detection of an ever- 
increasing number of high-z Quasi-Stellar Objects (QSOs), 
both radio-quiet and radio-loud. 

In Page et al. (2004a), the authors investigated the 
X-ray spectra of a sample of high-luminosity RQQs, find- 
ing that the continuum shape is essentially the same as for 
lower-luminosity (and, presumably, lower black-hole mass) 
radio-quiet AGN, i.e. Seyfert galaxies. The broad-band X- 
ray spectra could be well modelled either by a power-law 
(r ~ 1.9) together with one or two blackbody components 
to parametrize the soft excess, or by a double Compton- 



isation model (one component modelling the high-energy 
power- law; the other, the soft excess). Some of the RQQs 
showed evidence for neutral iron emission, although it has 
previously been found that the strength of the narrow Fe 
Kq line diminishes as the X-ray luminosity increases (Iwa- 
sawa & Taniguchi 1993; Nandra et al. 1997b; Page et al. 
2004b); this is the so-called X-ray Baldwin effect. In RLQs, 
these spectral features may be somewhat diluted, since some 
of the harder X-ray emission is thought to be Synchrotron 
or inverse Compton emission from the relativistic radio jet, 
rather than from the accretion disc. This also leads to RLQs 
tending to have flatter spectra over the 2-10 keV rest-frame 
(e.g., Williams et al. 1992; Lawson et al. 1997; Reeves et al. 
1997; Reeves & Turner 2000). 

This paper aims to analyse the X-ray spectra of a sam- 
ple of high-redshift QSOs, both radio-quiet and radio-loud, 
investigating, in particular, the continuum shape, iron emis- 
sion, Compton reflection and absorption. The results of such 
a comparison should help to determine just how fundamen- 
tally different RLQs and RQQs are: are the central engines, 
themselves, different, or are the differences solely due to the 
radio jet? A further major aim of the paper is to investi- 
gate how the X-ray spectral and continuum properties of 
high-redshift quasars (z = 2-6) compares to those at lower 
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Figure 1. A histogram showing the distribution of the redshifts 
for the RLQ, RQQs and BAL QSOs. 



redshifts. For instance, is there any evolution of the X-ray 
spectra with either luminosity or redshift, which might sug- 
gest that properties of quasar central engines evolve with 
time? Does the continuum photon index, or amount of in- 
trinsic absorption, vary, suggesting a changing environment? 

Section |21 explains how the data were processed and 
the spectra obtained, while Section |3 details the spectral 
analysis. The results are discussed in Section 0| 



carried out by Reeves et al. (2001). Before extracting the 
spectra, background light-curves were produced to identify 
any periods of high background. Most of the observations 
were found to show intermittent background flaring, so these 
time slots were subsequently ignored. However, the vast ma- 
jority of the observation of S5 0014-1-81 was affected by high 
background, so excluding all the flares was deemed to be im- 
practical. For this object, a very small region of 20 arcsec was 
used to extract the spectrum and, hence, minimise the back- 
ground contribution. Similarly small regions were used for 
PMN 0525-3343, PKS 2000-330 and RX J1028. 6-0844 in 
addition to excluding the worst periods of high background. 

For the other QSOs, spectra were then extracted from 
circular regions (40-50 arcsec in radius), with the back- 
ground spectra obtained from offset 'source-free' areas. Pat- 
terns 0-12 for MOS (0-4 for PN) were chosen, since there 
was no indication of pile-up. The SAS tasks rmfgen and AR- 
FGEN were then run, to produce the redistribution matrix 
functions and the ancillary response functions respectively. 
MOS 1 and MOS 2 spectra were co-added, after determin- 
ing there were no large discrepancies between them, and this 
combined spectrum was then fitted simultaneously with the 
PN data, using Xspec version 11.3, to produce the results 
given in this paper. 

Throughout this paper, a WMAP (Wilkin- 
son Microwave Anisotropy Probe) Cosmology of 
Ho = 70 km s"^ Mpc"\ SIa = 0.73 and Q.m = 
is assumed. Errors are given at the 90 per cent level (e.g., 
Ax^ = 2.7 for one degree of freedom). 



2 XMM-NEWTON OBSERVATIONS 

The sample in this paper constitutes all available XMM- 
Newton observations (until the end of 2004) of z>2 QSOs, 
with sufficient counts to allow spectral fitting (approxi- 
mately 500 counts in the combined MOS-I-PN results). 

Table 01ists the QSOs, their right ascensions and decli- 
nations, redshifts. Galactic absorbing columns (determined 
from the ftool nh; Dickey & Lockman 1990) and radio 
data, as well as information about the actual XMM-Newton 
observations. Where the NRAO/VLA Sky Survey (NVSS) 
did not detect the object at 1.4 GHz, the survey limit 
of 2.5 mjy is given, except for four cases (Q 0000—263, 
SDSS 1030+0524, SBS 1107+487 and APM 08279+5255) 
where more stringent limits were found in the literature. 
Figure plots the redshift distribution for the objects. 
The QSOs are a combination of proprietary data and pub- 
lic observations obtained from the XMM-Newton Science 
Archive^ . 

The ODFs (Observation Data Files) were obtained from 
the XMM Science Archive and processed using SAS v5.4^ for 
consistency, with the exception of PKS 0537—286. Because 
PKS 0537—286 was observed very early on in the mission 
(revolution 51), the ODFs are in a format which cannot 
be easily reprocessed. Hence, the results presented here are 
from the original data products; such an analysis was first 

^ http: / /xmm. vilspa.esa.es/external/xmm_data_acc/xsa/index.shtml 
^ The data for this paper were analysed during 2004, using the 
most recent calibration files available at the time. 



3 SPECTRAL ANALYSIS 
3.1 Simple power-law fits 

As is conventional, each of the spectra was first fitted with 
a simple power law, over the 2-10 keV (rest-frame) band; 
the resulting photon indices are given in Table |21 None of 
the QSOs showed statistically significant iron emission (the 
highest F-test probability was ~95 per cent), so only the 
upper limits for the equivalent widths of narrow lines are 
given in the table. For many of the sample including such a 
Gaussian led to unconstrained line energies; in these cases, 
the energy was fixed at 6.4 keV, corresponding to neutral 
iron. 

The highest redshift (z — 6.28) object, 
SDSS 1030+0524, showed a small amount of excess 
flux (one data bin in both the PN and MOS spectra) at 
approximately the correct energy for ionised iron emission, 
as previously noted by Farrah et al. (2004); Table |21 states 
the results for fitting neutral narrow lines only. Fitting the 
possible feature with a Gaussian gave an energy of 7.02'^^']^^ 
keV; this was approximately 92 per cent significant, with a 
high upper limit for the equivalent width of 3 keV. 

For each QSO, the rest-frame 2-10 keV power law fit 
was extrapolated over the full bandpass of XMM-Newton. 
Figure |21 shows the ratio plots formed by fitting a simple 
power-law over the observed frame 2-10 keV, and extrapo- 
lating over 0.3-10 keV for the RLQs, while Figures Inland 2] 
shows the results for the RQQs and BAL QSOs. From these 
three figures, it can be seen that the majority of the spec- 
tra showing curvature appear convex in shape, indicating 
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Figure 2. Power-law fits over the observed-frame 2-10 keV band are extrapolated over 0.3-10 keV, for the RLQs. In most cases this 
reveals evidence for curvature, generally in the sense of absorption at lower energies. MOS data are shown in black, PN, in grey. 
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Figure 3. As for Figure]^ but for the RQQs. 
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Figure 4. As for Figure |21 but for the BAL QSOs. Note the absorption edge just below 2 keV (observed) in the spectrum of 
APM 08279-f5255. 
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Table 1. Information about the QSOs and the XMM-Newton observations; the radio-loud and radio-quiet objects are separately ordered 
by Right Ascension. °- Galactic Nh, determined from the FTOOL nh; * from NVSS - NRAO/VLA Sky Survey (Condon et al. 1998). PKS 
0438-43 is outside the region covered; Rl = FsgHz/Fb; VignaU, Brandt & Schneider (2003); = Kuhn et al. (2001); ^ Turnshek et 
al. (1991); » Holt et al. (2004); (1) - values from Reeves & Turner (2000); (2) - 5 GHz from Gregory & Condon (1991); the B-band flux 
was calculated from NED; (3) - Elvis et al. (1994); (4) - Kaspi, Brandt & Schneider (2000); (5) - 5 GHz from Gregory et al. (1996); 
the B-band flux was calculated from NED; (6) — calculated from Zickgraf et al. (1997) and the NASA/IPAC Extragalactic Database 
(NED); (7) — for these objects, the 5 GHz upper limit was estimated from the 1.4 GHz value, using f^ oc u~^'^ for the radio power-law; 
the B-band flux was calculated from NED; (8) — calculated from Lonsdale, Barthel & Miley (1993) and NED; (9) - calculated from 
Parkes-MIT-NRAO (PMN) survey (Wright et al. 1996) and Fabian et al. (2001); (10) - 5 GHz upper limit estimated from the 1.4 GHz 
value; B-band flux from Bechtold et al. (2003); (11) - VignaU, Brandt & Schneider (2003); (12) - Kuhn et al. (2001); (13) - 5 GHz upper 
limit estimated from the 1.4 GHz value; B-band flux from Dai et al. (2004); (14) — 5 GHz upper limit estimated from the 1.4 GHz value; 
B-band flux from Grupe, Mathur & Elvis (2003) 



Object RA Dec Obs. ID Obs. date exposure time (ks) z Gal. Nj^^ 1.4 GHz log(RL)^ 

(J2000) (J2000) MOS 1 MOS 2 PN (lO^o radio 

cm-2) flux*" 
(mjy) 
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the possible presence of excess absorption. Some, however, 
are more concave, implying they may be better fitted by a 
broken power-law model. 

3.2 Complex continuum fits 

Many of the RLQs showed evidence for excess absorption 
in their rest-frames, between ~l-2 x IQ^^ cm~'^ (Table|2|l. 
Upper limits are given for the remaining objects. Using an 
ionised absorber model did not improve the fits further; 



however, the upper limits on the ionisation were not well 
constrained for most cases, so the presence of ionised ab- 
sorption cannot be excluded. This was even the case for 
APM 08279+5255 where an ionised absorption edge im- 
proved the fit. Figure|^plots the confidence contours for the 
absorbing column and photon index for the objects where 
excess Nh is detected. 

Broken power-law models were then fitted to the spec- 
tra of those QSOs with significant rest-frame absorption, 
fixing the Nh at the Galactic value. Also included were 
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Figure 5. 68, 90, 95 and 99 per cent confidence contours for the intrinsic absorbing column density and the photon index, F, for the 

QSOs where excess absorption is significant. 
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Table 2. Fits to the rest-frame 2-10 keV band. The line width was fixed at cr = 10 eV and, if the energy of the fitted line was 
unconstrained, it was fixed (f) at 6.4 keV; the EW is given for the rest-frame of each object. The F-test null probability is stated for the 
inclusion of the line; if there is no change in x^i then the F-test is undefined. Certain of the spectra were better fitted by the inclusion of 
excess (rest-frame) absorption; where this was the case, the value is given in the fifth column. " The spectrum of APM 08279-f 5255 also 

I Q 2g -1-0 22 

contains an ionised Fc K edge at 7.70_q'3j^ keV, with an optical depth, r = 0.44_q'2i. This object is also lensed, with a magnification 
factor of ~ 142, as are RX J0911-I-0551 (factor of 17) and HE 1104-1805 (factor of 12) - values from Dai et al. (2004); the luminosities 
given in this table (marked as **) are, however, uncorrected. The final column gives the limits on the reflection parameter; R = Q/2it, 
where Q is the solid angle subtended by the scattering medium. These were determined from fits to the full 0.3-10 keV XMM-Newton 
band. 
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PKS 0537-286 and 87GB 150844.6+571424, since the ra- 
tio plots in Figure |5| indicated the presence of curvature 
in these spectra. Using a broken power-law investigates 
whether what is being modelled by excess Nh is actually 
just the flattening of the continua at low energies, rather 
than real absorption. Tabled] gives the results of these fits 
and compares the reduced values with those found from 
the excess Nh fits. The final column shows which of the fits 
has the lower value. 

It was found that the spectra of S5 0014+81, RBS 
315, PKS 0438-43, PMN 0525-3343 and PKS 2126-15 
were all better fitted with excess absorption, while 
S5 0836+71, RX J1028.6-0844, RX J122135.6+280613, 
RX J0911.4+0551 and CSO 0755 were better modelled by 



a broken power-law. The two models gave the same re- 
suh for 87GB 142825.9+421804 and very similar values for 
RX J122135. 6+280613 and APM 08279+5255. Despite the 
statistical fit being slightly improved by using a broken 
power-law for these few objects, the values for Viow are very 
low and unphysical; hence, excess absorption appears to be 
a better explanation for these objects. 

Ahhough neither PKS 0537-286 nor 
87GB 150844.6+571424 required excess absorption, 
both spectra are noticeably better fitted with a broken, 
rather than a single, power-law. However, these continua 
curve in the opposite direction from the (possibly) absorbed 
sources, becoming steeper at lower energies. 

The existence of Compton reflection was also inves- 
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Table 3. Broad-band (0.3-10 keV, frame of observer) fits. The values given for Njj are in the rest-frames of the QSOs and are in 
additional to the Galactic values. Ax'^ gives the decrease in for the addition of the excess absorption. The final column gives the 
rest-frame energy bands corresponding to 0.3-10 keV in the observer's frame. The fit to APM 08279+5255 also includes an ionised 
absorption edge as detailed for the 2-10 keV fit. ** Again, the luminosities given in this table for APM 08279+5255, RX J0911+0551 
and HE 1104—1805 arc not corrected for the magnification by gravitational lensing (Dai ot al. 2004). 
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tigated in each QSO, after modelling the apparent ab- 
sorption at lower energies; the lack of iron line emis- 
sion suggests that no reflection components should be re- 
quired. Only PKS 0438—43 seemed to reveal evidence for 
a flattening of the spectrum at higher energies (Table I^J, 
with the other values being consistent with zero. Includ- 
ing the reflection component improved the fit statistic to 
X^/dof = 472/506 [F-test probability of 1.1 x IQ-^, com- 
pared to the power-law plus absorption model (Fit 4 in Ta- 
ble|21l]. The photon index then became P = 1.89 ± 0.06, with 
Nh = (1.46 ± 0.13) X 10^^ cm-2, for a refiection parameter 
of R = 0.62 ± 0.28. 

However, since the upper limit for iron emission in 
PKS 0438—43 is 33 eV, this precludes the presence of reflec- 
tion; hence, the deviation from the simple power-law must 
be due to intrinsic curvature of the spectrum. The broken 
power-law model was found to be an equally acceptable fit 
to the PKS 0438-43 spectrum (Table 



Summarising the work in this section, 14 of the 29 
QSOs can be fitted with a single power-law; of these, five 
are RLQs, seven RQQs and two BAL QSOs. The remaining 
15 objects show more complex continua, requiring either ex- 
cess absorption or a broken power-law. Nine of the 16 RLQs 
have intrinsic absorption above the Galactic value, with typ- 
ical columns around lO'^'^ cm-2. p^yj. q£ ^-^e six BAL QSOs 
are also absorbed (with X-ray columns of a few 10^^ cm-2, 
consistent with other BAL QSOs) while none of the RQQs 
shows evidence for excess absorption. 



4 DISCUSSION 

This paper presents XMM-Newton observations of a sample 
of 29 high-redshift (z>2), luminous (Lo.s-iokcV = 1 x 10''^ 
- 2 X lO**^ erg s-^) radio-loud and radio-quiet QSOs. Below 
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Table 4. Broken power-law fits to the spectra which show evidence for excess absorption; in these fits only Galactic N^f has been 
included. Fiow/rhigh signify the power-law slope fitting the low/high energy end of the spectrum respectively. The break energy is given 
in the rest frame of the object. The penultimate column restates the value obtained from a power-law plus excess rest-frame absorption 
model, first given in Table |21 except for the two marked *, for which no absorption was statistically required. The final column gives 
the change in between the fits. Here a positive value indicates that the x^ value is higher for the broken power-law compared to the 
absorption, meaning the fit is worse. 
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follows a discussion of the statistical properties of the high-z 
sample. 

4.1 2-10 keV 

Over the 2-10 keV rest-frame band, the mean power- 
law slopes were found to be: Frl — 1.55 ± 0.04; 
Frq = 1.98 to.ll; Tbal = 1.80 to.li- RLQs are known to 
have flatter slopes than the RQQs, with typical slopes of 
Trl ~ 1.6-1.7, compared with Frq ~ 1.9 (Williams et al. 
1992; Reeves et al. 1997; Reeves & Turner 2000; Page et al. 
2003; Piconcelh et al. 2005). The RLQs here have rather flat 
slopes, which could, conceivably, be an evolutionary effect. 
Bechtold et al. (2003) analysed a sample of high-z, radio- 
quiet AGN, finding that the average X-ray spectral index 
was flatter for the more distant objects, both due to the 
increased redshift and the higher luminosities. Grupe et al. 
(2005), however, found that T was steep (2.23 ± 0.48) for 
their high-z RQQs. Plotting 2-10 keV F against luminos- 
ity for the sample of objects in this paper agrees with the 
finding of Bechtold et al. (Figure |SJ , with weighted linear 
regression giving a slope of (—2.5 ± 0.3) x 10"'^. In this 
sample, however, the correlation only exists between F and 
luminosity, not redshift; this is consistent with there being 
no spectral evolution over time. 

This relationship between F and luminosity is, however, 
due to the dominance of the high- luminosity beamed, flat 
spectrum RLQs; considering only the RQQs and BAL QSOs 
in the sample, no correlation is present. The average photon 
index for the high-z RQQs in the current sample (1.98 to' 13) 
is very typical of low to medium-z AGN (Nandra et al. 1997a; 



George et al. 2000; Reeves & Turner 2000; Page et al. 2003; 
Porquet et al. 2004; PiconceUi et al. 2005) . This suggests (in 
agreement with earlier works, e.g.. Reeves & Turner 2000; 
Page et al. 2004a; Shemmer et al. 2005; Vignah et al. 2005) 
that there is no evolution of the spectrum with luminosity 
or redshift (for radio-quiet objects). 

An alternative method for obtaining a mean spectral 
shape is to fit all the data simultaneously, with a constant 
of normalisation to account for the flux difi'erences between 
the different spectra. Doing this leads to a photon index of 
F = 1.99 ± 0.05 (xVdof = 260/241) for the RQQ spec- 
tra (over 0.3-10 keV). Upper limits of EW < 108 eV and 
R < 0.53 for the equivalent width of a line at 6.4 keV and 
reflection parameter respectively were also derived. 

The photon index for the six BAL QSOs is formally 
consistent with that for the normal RQQs. Comparison with 
other BAL QSOs published in the literature (e.g., Wang et 
al. 1999; Gallagher et al. 2001; Green et al. 2001; Chartas, 
Brandt & Gallagher 2003; GaUagher et al. 2004) confirms 
that the mean slope found here is typical for such objects. 
This implies that the underlying X-ray continuum and cen- 
tral engine is the same as for RQQs, but the BAL QSO is 
being observed through the outfiow. 

Of the objects in this sample, none showed a significant 
Fe Ka line. The absence of iron lines is further confirma- 
tion of the X-ray Baldwin effect (Iwasawa & Taniguchi 1993; 
Nandra et al. 1997b; Page et al. 2004b), whereby the most 
luminous AGN, both radio-loud and radio-quiet, do not tend 
to reveal iron emission. In the case of the RQQs, many of 
the iron line equivalent widths are poorly constrained, be- 
cause such AGN are less X-ray bright (for a given optical 
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Figure 6. The more luminous QSOs tend to liave flatter 2-10 
keV (rest-frame) photon indices. Here the measured luminosi- 
ties for APM 08279+5255, RX J0911+0551 and HE 1104+0551 
have been decreased by factors of 142, 17 and 12 respectively to 
account for gravitational lensing effects (Dai ct al. 2004). The 
lower plot uses three luminosity bins (<1 X 10^*^, 1-10 X 10*®, 
>10 X lO'"^) to show this more clearly, with RLQs, RQQs and 
BALs all grouped together. 



luminosity) than their radio-loud equivalents, as mentioned 
in the introduction. The lack of reflection in the RLQs is 
not entirely surprising, since the radio-loud objects are be- 
lieved to be dominated by radio-jets. This means that the 
Doppler-boosted jet emission can be much stronger than any 
reflection component, thus diluting any features. 



4.2 Excess absorption 

13 of the 29 QSOs are better modelled by the inclusion of 
an additional absorption component; of these, nine (out of 
the 16 in the sample) are radio-loud, with four (out of six) 
being BALs. None of the seven radio-quiet objects require 
excess Nh. Radio-loud objects have often previously been 
found to be absorbed (e.g., Elvis et al. 1994; Brinkmann 
et al. 1997; Reeves & Turner 2000), while their radio- 
quiet counterparts are more seldom absorbed (Canizares 
& White 1989; Lawson & Turner 1997; Yuan et al. 1998; 
Page et al. 2003). Where Reeves & Turner (2000) anal- 
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Figure 7. Intrinsic column density (Njj) plotted against redshift. 
The Galactic absorption has been separately accounted for. RLQs 
are plotted as crosses/arrows, RQQs as circles/diamonds and 
BALs as stars/triangles, for detections/upper limits (90 per cent) 
respectively. 



ysed the same QSOs as those presented here, the val- 
ues found for Nh using the ASCA data were either in 
agreement with those from XMM-Newton (PKS 0438-43, 
S5 0836+71, 87GB 150844.6+571424, PKS 2000-330 and 
PKS 2126-15), or larger (S5 0014+81, PKS 0537-286, 
PKS 2149-306 and HE 1104-1805). RX J1028.6-0844 was 
also found to show strong (2 x 10^^ cm~^) absorption by 
ASCA (Yuan et al. 2000), while the XMM-Newton data pre- 
sented here, although absorbed, have a column density al- 
most a factor of ten lower. Also, a previous observation, 
reported in Grupe et al. (2004), find the spectrum to be 
largely unabsorbed. 

It is possible that this disagreement is due to the degra- 
dation of the ASCA SIS low energy response (see Yaqoob et 
al. 2000). The effect of this loss of response is an underesti- 
mate of the soft X-ray flux, which can manifest itself as an 
increase in Nh. However, since the ASCA and XMM-Newton 
observations are non-simultaneous, the absorbing columns 
could have changed. The data in Reeves & Turner (2000) 
show a mean column density of ~ 1.3 x lO'^^ cm~^, while 
the radio-loud quasars here give a very similar average of 
~ 1.4 x 10^^ cm~^ for the objects with detected Nh. 

Agreeing with earlier work using smaller samples of ob- 
jects (e.g., Elvis et al. 1994; Cappi et al. 1997; Reeves et al. 
1997; Fiore et al. 1998; Reeves & Turner 2000; Bassett et al. 
2004) found that the more distant QSOs were more heavily 
(intrinsically) absorbed - see their figure 9. To compare with 
the ASCA results, a similar plot was produced for objects 
observed by XMM-Newton. This includes absorbing columns 
for the QSOs in this paper, other datasets possessed by the 
authors and any other corresponding values obtained from 
the literature for objects at a redshift greater than 0.01. Fig- 
ure [7| shows the plot, while, in Appendix B, Table IHTl lists 
the objects, their redshift and the neutral absorbing column; 
for sources where there was evidence for a warm absorber, 
this was fitted before determining the additional cold col- 
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Because of the presence of censored (upper limit) data, 
the ASURV software (Astronomy Survival Statistics; Feigel- 
son & Nelson 1985; Isobe, Feigelson & Nelson 1986) was used 
to determine the line of best fit through the data points. Al- 
though this method takes into account that some of the mea- 
surements are upper limits, the errors on the actual values 
are not accounted for. The outcome of this analysis was that 
there is a positive correlation between redshift and intrin- 
sic Nh (linear regression slope of 0.61 ± 0.10); that is, the 
higher-redshift QSOs are more strongly absorbed. Kendall's 
Tau (also within ASURv) gave a probability of >99.9 per cent 
for the correlation. Spearman's Rank gave a probability of 
>99.99% for a positive correlation between the intrinsic Nh 
and redshift, both taking the upper limit Nh values as detec- 
tions and also setting them as the mid-points of the interval. 
The actual detections alone also showed this strong, positive 
correlation. 

The separate groups of RQQs and RLQs were also inves- 
tigated, to determine whether there is a difference between 
the intrinsic absorption in radio-quiet and radio-loud ob- 
jects. In Page et al. (2003) the authors compared a sample of 
XMM-Newton serendiptitous radio-quiet AGN with ASCA- 
detected radio-loud objects, concluding that there was a dif- 
ference between the absorption in these groups. The QSOs in 
that analysis were lower-redshift than those presented here, 
at z < 3. Using the two-sample tests within ASURV for the 
current sample gave a probability of 85-92 per cent for the 
RLQs and RQQs having different absorption, with Nh being 
higher for the RLQs. Although this difi'erence is marginal, it 
should be noted that none of the RQQs in this sample shows 
any evidence for an excess column, while more than half of 
the RLQs are absorbed. These results are not conclusive, 
though they do tend to indicate that there is a difference 
in the intrinsic absorbing columns found in radio-loud and 
radio-quiet AGN. Recently, however, Grupe et al. (2005) 
found that their high-redshift RLQs were not more highly 
absorbed than the radio-quiet counterparts; only some of 
their targets overlap with the sample prsented here. Clearly 
work remains to be done in this area, with more good quality 
observations of high-z quasars, both radio-loud and radio- 
quiet, required. 

It should be noted that, since RLQs can, in general, 
be studied out to higher redshifts (see the introduction), 
the Nh-z correlation and the fact that RLQs may be more 
strongly absorbed than RQQs could be due to the same ef- 
fect. The Partial Spearman Rank test can be used to try and 
determine whether, for example, the correlation between Nh 
and Rl is simply due to high-redshift quasars frequently be- 
ing radio-louder than the nearby ones; in this case, the 'real' 
correlation would be Nh-z. Unfortunately, because many of 
the estimates for the absorbing column are upper limits, as 
are the Rl values for the RQQs, it was not possible to ob- 
tain a definitive answer on this subject. Taking the upper 
limits as actual measurements gave similar probabilities for 
the Nh-z and Nh-Rl correlations, while halving each of the 
upper limits gave a slightly higher probability that Nh is 
correlated with redshift only due to its correlation with Rl. 
In each case, however, individual Spearman Rank tests for 
Nh-z and Nh-Rl implied that Nh was more strongly corre- 
lated with redshift. 



Table 5. The QSOs for which intervening Damped Lyman-Q 
systems have been firmly identified or discounted. The reference 
for the information is given in the final column. 



Object 


DLA ? 


XMM Nh 
(10^2 cm-2) 


DLA reference 


S5 0014-1-81 


no 


1.82 ± 0.19 


Bechtold et al. 








(1994) 


Q 0420-388 


yes 


<0.76 


Elvis et al. (1994) 


PMN 0525-3343 


no 


1.91 ± 0.38 


Peroux et al. 








(2001) 


PKS 0537-286 


yes 


<0.04 


Kanekar & Chen- 








galur (2003) 


RX J1028.6-0844 


yes 


<1.17 


Peroux et al. 








(2001) 


Q 0000-263 


yes 


<0.82 


Lu et al. (1996) 


LBQS 0053-0134 


no 


<1.54 


Storrie- Lombard! 








& Wolfe (2000) 


HE 1104-1805 


yes 


<0.13 


Lanfranchi & 








Friaga (2003) 


BR 2237-0607 


yes 


<3.77 


Lu et al. (1996) 



4-2.1 Damped Lyman-a systems 

It is possible that some of the neutral absorption identified in 
these spectra is not intrinsic to the QSOs, but, rather, is due 
to intervening damped Lyman-a (DLA) systems. DLAs are 
high density (Nhi > 2 x lO'^^ atoms cm~^) systems, which 
contain most of the neutral gas found at high redshifts; the 
Lyman-a forest lines are far more numerous, but only make 
up a small fraction of the gas. Searching the literature, it 
was found that some of the QSOs in the current sample do 
have intervening DLAs. However, these do not tend to be 
the objects which show the highest excess absorption (see 
Table 1^. Hence, most of the absorption in these spectra is 
likely to be intrinsic to the QSOs themselves. 

4-2.2 Origins for the absorbing column 

Assuming that the absorption is intrinsic to the objects, 
there are various possible origins for the absorbing material 
(see, e.g., Elvis et al. 1994; Cappi et al. 1997). The quasars 
themselves could be surrounded by a region of dust or gas 
(Rieke, Lebofsky & Kinman 1979; Sanders et al. 1989; Web- 
ster et al. 1995), which would lead to obscuration of the 
AGN. Alternatively, since some RLQs can be located in 
clusters of galaxies (Yee & EUingson 1993), the absorbing 
material measured may be due to a cooling fiow (White et 
al. 1991; Elvis et al. 1994). Cold clouds of material close 
into the central regions of the AGN may exist (Guilbert & 
Rees 1988; Celotti, Fabian & Rees 1992), which could lead 
to both reflection features and absorption. That being said, 
no Compton reflection humps or iron lines are found in the 
present sample. 

None of the absorbed (non-BAL) quasars have signifi- 
cantly reddened optical spectra, indicating that the dust/gas 
ratio is low. The absorbing material should, therefore, be lo- 
cated close to the central engine, in order for the dust to 
sublimate. 

The possibility exists that the material is related to the 
warm absorbers seen in Seyfert Is; the column densities of 
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~ 10^^ cm~^ are similar to a number of these, although 
the level of ionisation cannot be tied down. Because of the 
bandpass covered for these high-z quasars, it may be just the 
continuum curvature below the iron K band which is being 
observed. Examples of this (convex) spectral curvature arc 
seen in NGC 3783 (Kaspi et al. 2002; Reeves et al. 2004) and 
NGC 3516 (Turner et al. 2005). In order to constrain such an 
absorber, lines and edges would need to be measured at high 
resolution and signal-to-noise, which will not be possible for 
these sources until the next generation of instruments (e.g., 
Con-X and XEUS). 

If RLQs do tend to be more absorbed than RQQs, then 
it is possible that the absorbing material is somehow linked 
to the radio jet - perhaps some kind of confining medium, 
which helps to coUimato the jets (e.g., Crawford & Vander- 
riest 1997; Ferrari 1998). 



5 SUMMARY 

As has been found in previous samples (Williams et al. 1992; 
Reeves et al. 1997; Reeves & Turner 2000; Page et al. 2003; 
Piconcelli et al. 2005), RLQs have flatter spectral indices 
than RQQs, suggesting the presence of an extra component 
due to a relativistic jet. The lack of any iron emission lines or 
a reflection component in the RLQs is also consistent with 
this interpretation. However none of the radio-quiet quasars 
shows significant Fc emission either, but this is probably due 
to the poorer photon statistics in these objects. 

There is no evidence for any evolution of the X-ray 
continuum for radio-quiet quasars with either redshift or 
luminosity. The mean photon index for this high z, high 
luminosity sample is consistent with those in lower z or 
lower luminosity samples (Nandra et al. 1997a; George et 
al. 2000; Reeves & Turner 2000; Page et al. 2003; Porquet 
et al. 2004; Piconcelli et al. 2005). Thus the X-ray proper- 
ties of the quasars do not appear to evolve with redshift or 
luminosity and may just simply scale with black hole mass 
and luminosity. 

The continuum properties of six high-z BAL QSOs have 
also been measured in this paper. Four of these six ob- 
jects show large absorbing columns (>10^^ cm~^ in the rest 
frame), presumably due to intersecting the outflowing mar- 
terial along the line of sight. After correcting for absorption, 
the continuum photon indices of the BAL QSOs are the same 
as the non-BAL RQQs, consistent with the difference solely 
being due to the inclination of the observer with respect to 
the outflow. 

Over half of the RLQs in this sample (nine out of 16) 
show evidence for intrinsic absorption in excess of the Galac- 
tic value. In the QSO rest frame the absorption is typically 
10^^ cm~^, in agreement with earlier studies with ROSAT, 
ASCA and Chandra (Elvis et al. 1994; Brinkmann et al. 
1997; Reeves & Turner 2000; Bassctt et al. 2004). There is 
no correlation between the absorption and the occurence of 
any damp led Lya system along the line of sight, suggest- 
ing that the absorbers are associated with the quasar or its 
host galaxy. The lack of reddening in the optical spectra of 
the quasars favours the absorber being located close (within 
a few parsec) of the quasar central engine. Absorption is 
not found in low-z radio- loud quasars [e.g., 3C 273 (Page et 
al. 2004d); B2 1028-1-31, B2 1128-1-31, B2 1721-h34 (Page et 



al. 2004c)], consistent with the N/f-z correlation found from 
ASCA quasar samples (e.g., Reeves & Turner 2000). High 
throughput and high resolution X-ray spectroscopy will be 
needed to determine the origin of the low energy absorp- 
tion, through the detection of absorption lines/edges in the 
quasar spectra. 
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APPENDIX A: PREVIOUS WORK 

The data for some of these XMM-Newton observations 
has been previously published. As already mentioned, 
PKS 0537-286 was analysed in Reeves et al. (2001), with 
the current results being in agreement. 

Worsley et al. (2004a) published a paper about the ab- 
sorption in PMN 0525—3343, using data from eight separate 
observations. Their value for the intrinsic Nh (for a joint fit 
to all the datasets, rather than the single observation consid- 
ered in this paper) is consistent with the results here at the 
90 per cent level. Although there are no clear indications of a 
warm absorber, they conclude that cold absorption is incon- 
sistent with the optical data and that an ionised absorber 
would be more likely. Both of the XMM-Newton observations 
for 87GB 142825.9-1-421804 were also analysed by Worsley 
et al. (2004b), while only the second (longer) dataset is used 
here. Once again, they favour a warm absorber model, al- 
though intrinsic cold absorption is not ruled out. The column 
density for cold absorption is in agreement with the value 
presented in this paper. 
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Ferrero & Brinkmann (2003) analysed the spec- 
tra of PKS 2126-15, PKS 2149-306, Q 0000-263 
and Q 1442+2931. The evidence for absorption in 
PKS 2149—306 is marginal and only the Galactic column 
is required for Q 0000-263 and Q 1442+2931; these results 
are in agreement with the current work. PKS 2126—15 is 
found to be absorbed, with neutral Nh ~ 1.4 x 10^^ cm~^; 
this is slightly larger than the value found here, but their 
spectrum was fitted down to 0.2 keV, rather than 0.3 keV. 

Brocksopp et al. (2004) presented the spectra of 
PKS 0438-436, RX J122135.6+280613, PKS 2000-330 and 
SBS 1107+487, finding that the first two required cold 
absorption components (with values consistent with those 
given here), while the other two were not significantly ab- 
sorbed (also found here). 

Yuan et al. (2004) analysed the observation of 
RX J1028.6— 0844 discussed in the present paper. An ear- 
lier observation (Obs. ID; 0093160701; date of observation: 
2002-05-15) was covered as well and is also the subject of 
the paper by Grupe et al. (2004). This previous observation 
was much shorter (<5 ks in each of the EPIC instruments), 
so the data are not included in this paper directly; how- 
ever, they were analysed and the results were in agreement 
with Grupe et al. (2004), in that no intrinsic absorption was 
found. In the more recent observation, however, there is ev- 
idence for excess Nh, as also found by Yuan et al. (2004). 
They considered the MOS and PN cameras separately, find- 
ing that the values of Nh from the PN detector were lower 
than those from MOS - possibly underestimated due to low- 
energy calibration uncertainties. 

Hasinger, Schartel & Komossa (2002) discovered the 
same ionised Fe K edge in the XMM-Newton spectrum of 
APM 08279+5255 as was found in this present analysis. 



Table Bl. The data for Figure |7| listing the corresponding red- 
shifts and neutral Nh values; the objects are listed by decreasing 
redshift. * marks the objects which also have a warm absorber. 
All values, except where the quasars are marked were inferred 
from the authors' own work. Values from the literature: " Por- 
quet et al. (2004); " Gallo et al. (2004b); Gallo et al. (2004a); " 
Vaughan et al. (2004). 
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HE 1104—1805 


2.305 


<0.13 


LBQS 0053-0134 


2.062 


<1.54 


PG 1247+267 


2.038 


<0.09 


PB 5062 


1.77 


0.52 ± 0.32 


SBS 0909+532 


1.376 


<0.01 


PG 1634+706 


1.334 


<0.07 


PHL 1092 


0.396 


<0.001 


Mrk 493 


0.313 


<0.005 


UM 269 


0.308 


<0.08 


Q 0144-3938* 


0.244 


<0.212 


PG 0953+414t'^ 


0.2389 


<0.006 


PG 1427+480^'^ 


0.221 


<0.016 


PG 0947+396^'^ 


0.206 


<0.007 


PG 1116+215^''^ 


0.177 


<0.019 


PG 1322+659^'^ 


0.168 


<0.028 


PG 1048+342 


0.167 


<0.03 


PG 1202+281^'^ 


0.165 


<0.005 


PG 1402+261^'^ 


0.164 


<0.028 


Mrk 1014 


0.163 


<0.02 


PG 1307+085^'^ 


0.155 


<0.07 


PG 1115+407^'^ 


0.154 


<0.008 
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Table Bl. - continued 



Object 


redshift 


AT /1 n22 — 2\ 

Nh (10 cm ■') 


PG 1352+183+'" 


0.152 


<0.007 


PG 1114+445* 


0.144 


0.13 ± 0.06 


PG 1626+554 


0.133 


<0.02 


Mrk 876 


0.129 


<0.07 


IH 0419-577 


0.104 


<0.04 


PG 0804+761 


0.100 


<0.06 


HE 1029-1401 


0.086 


<0.008 


Mrk 1383 


0.086 


<0.01 


PG 1211+143 


0.0809 


<0.004 


Mrk 478t''' 


0.077 


0.022 ± 0.011 


Mrk 205 


0.071 


<0.05 


Mrk 304* 


0.066 


<0.05 


MR 2251-178 


0.064 


0.08 ± 0.02 


PG 0844+349 


0.064 


0.03 ± 0.02 


IZwlt." 


0.0611 


0.05 ± 0.006 


PG 1244+026t''' 


0.048 


0.121 ± 0.001 


Fairall 9 


0.047 


0.08 ± 0.05 


MCG -2-58-22 


0.047 


<0.05 


IH 0707-495t'<* 


0.0411 


0.16 ± 0.002 


ESQ 141-G55 


0.036 


<0.02 


Mrk 841 


0.036 


<0.001 


Mrk 509 


0.0344 


0.011 ± 0.004 


Ark 120 t'*^ 


0.0323 


<0.02 


Mrk 335 


0.026 


<0.02 


Mrk 896 


0.026 


<0.01 


iVirK oDy 


U.Ul 1 




Mrk 1044 


0.017 


<0.005 


NGC 5548 


0.017 


<0.001 


IC 4329a* 


0.016 


0.16 ± 0.04 


Mrk 766* 


0.0129 


0.011 ± 0.003 


BAL 






Q 0000-263 


4.098 


<0.82 


APM 08279+5255 


3.87 




CSO 0755 


2.86 


1.52 ± 0.38 


RX J0911.4+0551 


2.80 




SBS 1542+541 


2.371 


<1.16 


Q 1246-057 


2.236 


1.0D_Q g2 



